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Abstract. We report measurements of the 12 C/ 13 C abundance ratio in the three galactic regions G 333.0-0.4, 
NGC 6334 A and G 351.6-1.3 from observations of the 12 Cl 3 P2 ^ 3 Pi transition and the hyperfme components 
of the corresponding 13 Cl transition near 809 GHz. These transitions were observed simultaneously with the CO 
7-6 line emission at 806 GHz with the AST/RO telescope located at the South Pole. From a simultaneous fit to 
the 12 Cl 3 P2 —* 3 Pi transition and the HF components of the corresponding 13 Cl transition and an independent 
estimate of an upper limit to the optical depth of the 12 Cl emission we determine intrinsic 12 Cl/ 13 Cl column 
density ratios of 23 ± 1 for G 333.0-0.4, 56 ± 14 for NGC 6334 A and 69 ± 12 for G 351.6-1.3. As the regions 
observed are photon dominated, we argue that the apparent enhancement in the abundance of C towards G 
333.0-0.4 may be due to strong isotope-selective photodissociation of 13 CO, outweighing the effects of chemical 
isotopic fractionation as suggested by models of PDRs. Towards NGC 6334 A and G 351.6-1.3 these effects appear 
to be balanced, similar to the situation for the Orion Bar region observed by Keene et al. (1998). 
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1. Introduction 

The study of CNO isotope abundance ratios in the in- 
terstellar medium is crucial to understanding galactic 
chemical evolution. In particular, the ubiquity of C-based 



molcciiles in interstellar clouds has made the 1J C/ ld C 
abundance ratio an important chemical diagnostic. From 

measurements of C 18 and 13 C 18 in 13 interstellar 



clouds, Langer fc Penzias (199C, 1993) found a galactic 



gradient in the local carbon isotope abundance ratio (here- 
after: n(C)-ratio) ranging from 25 towards the galactic 
center to about 60 — 70 in the local ISM and out to a 
galactic radius of about 10 kpc. They also found an in- 
crease in the n(C)-ratio derived from CO towards clouds 
exposed to higher UV radiation fields (namely Orion KL 
and W33), supporting the suggestion made by models of 
photon dominated regions (PDRs) that depletion in 13 CO 
can be ascribed to inefficient self- shielding resulting in 
isoto pe 
19881 ) 



As discussed by Keene et al. (1998 ), who report the 
first detection of the strongest of the 13 Ci hyperfme (HF) 
structure components at 809 GHz from the Orion Bar 
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region, the abundance of atomic carbon is sensitive to 
the effects of chemical isotopic fractionation and isotope- 
selective photodissociation, allowing the study of the sig- 
nificance of these competing effects in PDRs and to verify 
model predictions (e.g. Lc Bourlot et al. 1993| , Koster et al 



1994). The charge exchange fractionation reaction, 13 C 4 
T^O ^ 12 C+ + 13 CO + 36 K, being exothermic, pref- 
erentially incorporates 13 C + into 13 CO, thus leading to a 
local 13 C + and 13 G depletion and a corresponding 13 CO 
enhancement. This may be balanced by isotope-selective 
photodissociation, which in contrast reduces the gas-phase 
abundance of the less efficiently self-shielded 13 CO iso- 
topomer ( van Dishoeck fc Black 19"8§| ) . Models show that 
the proportion of these effects strongly depends on the 
temperature derived for the C+/CI/CO transition zone 
(due to the low energy barrier of the exchange reaction) 
and the comparison between different models gives non- 



selective photodissociation jvan Dishoeck fc Black P993) 



conclusive results (Koster et al. 1994; Le Bourlot et al 



From their measurement of the 12 G 3 P2 — > 3 Pi and the 
strongest HF comp onent of the 13 Ci equivalent transition, 
|Keenc et al. (1998|) find a N( 12 C)/N( 13 C) column density 
ratio (hereafter: N(C)-ratio) of 58 ± 12 towards a position 
near the western end of the Orion bar. They derive a some- 
what higher ratio of 75 ± 9 towards the same position from 
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observations of C 18 and 13 C 18 0. In comparison, direct 



observations of the 12 C + / 13 C + abundance ratio (Stacey 
et al 



993 



Borciko & Betz 1996), which should not be 



significantly influenced by chemical isotopic fractionation 
or isotope-selective photodissociation, yield a n(C)-ratio of 
58±6. Keene et al. (1998) conclude that, in contrast to the 
model results by Lc Bourlot et al. (1993) the importance 
of chemical isotopic fractionation is (almost) compensated 
for by the isotope-selective photodissociation. Whether 
this compensation is peculiar to the special conditions of 
the Orion Bar, or whether it also holds for other massive 
star-forming cores is a significant question when address- 
ing galactic chemical evolution, since our current under- 
standing is almost entirely based on the observed carbon 
monoxide isotopomer abundance ratios. It is therefore im- 
portant to extend the study of 13 Cl to other regions. 

In this Letter, we report measurements of the 
3p 2 _>3 Pi trans ition of 12 Ci at 809341.97 MHz and the 
F= 5/2 - 3/2, F= 3/2 - 1/2, and F= 3/2 - 3/2 com- 
ponents of the equivalent transition of 13 Cl at 809493.7 
MHz, 809125.5 MHz, and 809121.3 MHz, respectively, 
towards three galactic star-froming regions associated 
with strong IRAS sources and CO emission: G 333.0- 
0.4, part of the large H n complex RCW 106 and as- 
sociated with a molecular cloud with bright CO lines 
( |Gillcspie et al. 1977 , Brand et al. 1984) at a distance 



cc ntinuum region (McBreen et al. 1973) associated 
bright PDR ( [Burton et al. 2000] , 



1 seated within a giant molecular cloud at a distance 



of 30 i 



Kraemer et al 



of 4.2 frpc flShaver fc Goss 197C | ); NGC 6334 A, a bright 
FIR 
with a 
200C()J 
of 1.7 

ded O 
5 kpc 
selected 
sition 



kpc (Neckel 1978); and G 351.6-1.3, a luminous 
1O 6 L ) co mpact H H region exci ted by an embed- 
B cluster (McBreen et al. 1983 ) at a distance of 
( Radhakrishnan et al. 1973| ). These regions were 



because of strong e mission in the 3 Pi ^P tran- 



of 12 Cl at 492 GHz QHuang et al. 1999|) from a list 



mdidatcs of galactic II II regions having st rong C O 



(Brand et al. 1984, Hcnning fc Launhardt 1998) and CS 



( Bronfman, Nyman, fc May 1996 , Launhardt et al. 1998 ) 
emission lines. 
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Broad band spectra (50 hrs of integration time 
e) towards the galactic star-forming regions NGC 
;top), G 333.0-0.4 (middle) and G 351.6-1.3 (bot- 
e insets show the same spectra magnified in y- 
le 12 Cl and 13 Cl lines were observed in the lower 
simultaneously with the CO 7-6 line in the up- 
side|band. Note the erratic baseline toward the 13 Cl 
in G 351.6-1.3, which resulted in a skewed 
fit and precludes a clear 3cr assignment of this 



Th 



enn ssion 



(1998) (cf. Cologne Database for Molecular Spectroscopy 
- www.cdms.de) were used to determine the line centroids. 



2. Observations 



3. Results 



The observation were made in the Austral win- 
ter of 2000 with the 800 GHz receiver at at 
the Antarctic Submillimctcr Telescope and Remote 
Observatory (AST/RO, cf. |Stark et al. 200l| ). At these 
frequencies the FWHM beam size is « 80". As back- 
end we used an AOS with a velocity resolution of 0.25 
km s _1 at the observed frequencies. The zenith opacities 
at 809 GHz at the times of observation were « 1.0 ± 0.2 
System temperatures on average were between 11,000 K 
for G 333.0-0.4, 18,000 K for NGC 6334 A, and 25,000 
K for G 351.6-1.4. The 12 Ci and 13 Ci lines were observed 
in the lower sideband simultaneously with the CO 7-6 
l ine in the upper s ideband. The spectra were calibrated 
( jBtark et al. 2001 ) to the effective radiati on tempera- 
ture Tr* scale. The frequencies reported by Klein et al 



Fig. 1 shows the spectra obtained towards the three re- 
gions. The strongest HF component of 13 C at 809 GHz, 
F= 5/2 - 3/2, is clearly visible in the G 333.0-0.4 and 
NGC 6334 A spectra; it is marginally detected in G 351.6- 
1.3. The weaker HF structure satellites, F=3/2 — 1/2 and 
F=3/2 — 3/2, are still hidden in the noise in all cases. 



Following the first detection towards Orion by Keene et al 
(1998j ), these regions are the only other 13 C detections to 



date. We derive the 12 Cl/ 13 Ci intensity ratio by a simul- 
taneous Gaussian fit of the lines (including all 13 Cl HF 
components) with fixed spacing and a single line width, 
and with the ratio of the 13 Cl amplitudes fixed to their 
quantum mechanical values of 0.600:0.333:0.067; the free 
fit parameters are the common width and LSR-velocity, 
a common amplitude and the relative line intensity ratio 
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Table 1. Line fit results 



5 Pi -^ 3 P emission lines. The 492 GHz 12 Ci 3 Pi -^ 3 P 



emission 
line 



Tr* 
[K] 



Av ls 



[kr 



a 

ha 



NGC 6334 A 



CO 
CO 



0.10 ±0.02 4.9 ±0.2 
8.04 ±0.06 5.3 ±0.1 
0.14 ±0.03 32 ±7 



45 ± 11 



G 333.0-0.4 



^CT 

13 Ci 
CO 
CO 



0.12 ±0.02 
0.04 ±0.01 
3.57 ±0.02 



5.0 ±0.1 

18 ± 1 
5.4 ±0.1 



18 ±1 
45 ±4 



narrow 
broad 
n 



0.46 ±0.02 12.6 ±0.4 



G 351.6-1.3 



-is 



Ci 



0.05 ±0.01 4.4 ±0.1 55 ±10 



CO 
CO 



2.15 ±0.02 
1.77 ±0.02 



3.3 ±0.1 
8.9 ±0.2 



line transition has also been measured (Huang ct al. 1999) 
towards these sources; assuming that the emission is ex- 
tended and need not be corrected for different beam sizes 
at 492 and 809 GHz, we derive a line intensity ratio of 
0.48 for G 333.0-0.4, 0.63 for NGC 6334, and 0.75 for G 
351.6-1.4. If the sources are centrally peaked, coupling to 
the narrower 809 GHz beam would be higher, and the 
intrinsic line ratios would be correspondingly lower. 

We use a single component escape probability excita- 



tion code (Stutzki & Winnewisser 1985) to interpret the 



The 13 Ci amplitudes given are the sum over the three HF 
components (see text). 



observed line intensities. As an additional constraint we 
need an estimate of the sour ce density. Only for NGC 
6334 do we have an estimate: Kraemer, Jackson, & Lane 
(1998|) compared their observed line intensities of C n 
(158/im) and Oi (146/im & 63/im) and the CO line and 



FIR continuum intensities to PDR models of Wolfire 



Tielens, fc Hollenbach (1990 ) to estimate the gas den- 



a = fy 3 -. We included the CO 7-6 line from the other 
sideband as an additional, independent Gaussian. Due 
to broad line wings apparent in the CO and 12 Ci emis- 
sion lines towards G 333.0-0.4, we fit two line components 
(broad & narrow) for all emission lines here. In the case 
of NGC 6334 A and G 351.6-1.4, where no wings are ap- 
parent for the 12 Ci line but asymmetric line profiles are 
conspicuous for the CO 7-6 line, we fit a second broad line 
component to the CO emission line only. The fit results for 
all sources are given in Table 1. The uncertainties quoted 
are the formal fit errors based on a la excursion. Since 
the 13 Ci and 12 Ci lines were measured simultaneously in 
the same receiver sideband, calibration uncertainties are 
negligible. 

In the optically thin case, the N(C)-ratio is directly 
given by the line intensity ratio a; beyond this limit an 
optical depth correction must be applied. We prefer this 
approach to the alternative of including the optical depth 
into the fit by fitting saturated line profiles of the form 
1 — cxp ( — t x </>„), because of the macro-turbulent nature of 
the velocity field in the ISM and the corresponding failure 
to observe saturated profiles even for optically thick lines 
of 13 CO in any astronomical source. This is supported by 
the Gaussian profile fit to the CO 7-6 line with a width 
only marginally wider than that of 12 Ci. 



3.1. Optical Depth Estimate 

To derive the intrinsic 12 Ci/ 13 Ci column density ratio, 
N(C), we apply an optical depth correction to a. Following 
Keene et al. (1998| ), the ratio is given by 



A^( 12 C?) _ a 



(1) 



with the line escape probability /3(t) = (1 — e~ T )/r. To 
further constrain the properties of the Ci emitting gas, 
we use the line intensity ratio of the 12 Ci 3 P 2 ^ 3 Pi and 



sity and FUV field towards the embedded radio contin- 
uum sources. Towards NGC 6334 A they derive an aver- 
age FUV field of log Go = 4.5 and a total gas density of 
n(H2) = 6.3 x 10 3 cm~ 3 within their 1' beam. It is fair 
to assume that the density of the Ci emitting gas in NGC 
6334 A and the other two sources is similar or higher, as 
it is located deeper in the molecular cloud. Densities of 
that magnitude are also supported by the bright CO 7-6 
emissions we observe. 

As is shown in Fig. 2, this density immediately con- 
strains the temperature of the Ci emitting region to 
about 40 K or lower: at higher temperatures, the 12 Ci 
3 P 2 — > 3 Pi/ 3 Pi — > 3 Po line ratio at these densities would 
be substantially higher than observed. This temperature 
estimate is perfectly consistent with the scenario of ex- 
tended PDR emission. In this regime of density and tem- 
perature, the observed 12 Ci 3 P 2 -^ 3 Pi/ 3 Pi -^ 3 Po line ra- 
tio constrains the column density per unit velocity interval 
to values below 2 x 10 17 cm -2 / km s _1 and corresponding 
optical depths of the 12 Ci 3 P 2 — > 3 Pi line below r w 0.4. 
At higher column densities trapping would push the line 
ratio to much higher values than observed. 

The total column density in the « 5 km s _1 wide 12 Ci 
lines is then about 10 18 cm -2 ; this corresponds, with the 
standard gas phase abundance of carbon of about 10~ 4 , to 
a hydrogen column density of 10 22 cm~ 2 . Considering that 
each PDR surface has a Ci column density corresponding 
to an Ay of a few, we therefore conclude that we see a few 
PDR surfaces per beam, consistent with the usual scenario 
of a clumpy, UV-penetrated massive cloud core. 

4. Discussion 

For the range of density, temperature and column density 
derived above, our radiative transfer model predicts line 
brightnesses for the 12 Ci 3 P2 ^ 3 Pi lines of around 5—10 
K. Since the measured line brightnesses are a few degrees 
in the observed sources, we estimate a beam filling factor 
close to unity, consistent with extended, smooth emission 
within the AST/RO beam. This scenario suggests the 12 Ci 
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in the PDR models of Lc Bourlot et al. (1993), which 



t = 0.4 



ot 30 K 
at 40 K 
at 60 K 



2 3 4 5 log n(H 2 ) 6 

Fig. 2. Results from the line escape probability code: Bold 
lines show 12 Cl line ratios for 40 K. At these temperatures 
and lower (30 K, dashed), densities of n(H 2 ) = 6.3 x 10 3 
cm~ 3 indicate an optical depth of r < 0.4 and constrain 
the densities to N( 12 C) < 2xl0 17 cm~ 2 /kms~ 1 (indicated 
by a cross). At higher temperatures (60 K, dotted) the 
ratios at these densities would be substantially higher than 
observed. 



emission is moderately optically thin. For the rough up- 
per limit of the optical depth r w 0.4, the optical depth 
correction factor l//3(r) = 1.21 only raises the N(C)-ratio 
derived above by at most 25%. We estimate the N(C)-ratio 
in G 333.0-0.4 to be 23 ± 1 or lower (56 ± 5 for the weak 
extended component), that in NGC 6334 A to be 56 ± 14 
or lower, and that in G 351.6-0.4 to be 69 ± 12 or lower. 

The N(C)-ratios towards G 351.6-1.4 and NGC 6334 
A are consistent with the intrinsic isotopic ratios and do 
not call for significant isotopic fractionation or isotope- 
selective photodissociation. These results confirm the con- 
clusion of Keene et al. (1998| ) for their Orion Bar data. 
Towards G 333.0-0.4, however, we find the ratio to be 
substantially lower than the average values derived from 
CO isotopomers. In this context one should recall, that 
(Langer & Penzias 1993) find an increase in the CO/ 13 CO- 
isotopomer ratio towards clouds exposed to higher UV ra- 
diation fields, supporting the suggestion of PDR models 
that the depletion in 13 CO can be ascribed to less effective 
self-shielding and enhanced isotope-selective photodisso- 
ciation ( van Dishoeck fc Black 1988| ). This would then 
yield an increase in the abundance of the 13 C isotope, as 
is shown in the PDR-models by Koster et al. (1994), if 



not counter-balanced by chemical isotopic fractionation. 
As chemical isotopic fractionation becomes more efficient 
with decreasing temperatures and eventually dominates 
over isotope-selective photodissociation, there will be a de- 
crease in the abundance of the 13 C isotope. This is shown 



produce a lower temperature in the transition zone due 
to reduced heating. Cloud surface temperatures therefore 
influence the balance of the two competing effects and, 
ultimately, determine the carbon isotope abundance ra- 
tio. The temperature structure in PDR surfaces, details 
of which are not consistently explained by present PDR 
models, might be responsible for the very different atomic 
carbon isotope abuncance ratios observed to date. 
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